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1 .  INTRODUCTION 


The  naterial  Gd^SCjCa^O^  is  a  mixed  garnet  with  the  same  crystal  struc¬ 
ture  as  the  popular  laser  material  Y3A15C>12  (YAG).1'2  The  chemical  formula 
for  YAG  is  sometimes  written  as  Y3A12A13C>12  to  illustrate  that  two  types  of 
aluminum  sites  exist  in  this  material.  In  the  material  Gd3Sc2Ga3012  each  ion 
ty^e  occupies  a  single  site  type.  The  optical  spectra  and  laser  action  of 
Nd  in  Gd3Sc2Ga30^2  have  been  reported  by  Kaminskii.3'4  An  interesting 
feature  of  this  material  is  that  if  it  is  doped  with  both  Nd3  +  and  Cr3+, 
energy  is  transferred  efficiently  from  the  Cr3  +  to  the  Nd3  +  impurities.5  The 
microwave  phonon  attenuation  of  this  and  related  garnets  has  been  measured  at 
a  number  of  different  temperatures.5  The  indices  of  refraction  of  a  number  of 
related  mixed  rare-earth  garnets7  have  been  measured  in  the  wavelength  region 
from  5200  to  6200  A.  Although  the  optical  spectra  of  only  one  rare-earth  ion 
in  this  host,  Nd3+,  have  been  reported,  we  considered  the  host  important 
enough  to  do  a  preliminary  crystal-field  analysis  of  rare-earth  ions. 

In  this  report  we  use  techniques  reported  previously8'9  to  obtain  first 
estimates  of  the  crystal-field  parameters  for  Nd: Gd3Sc2Ga301 2»  These  parame¬ 
ters  are  used  as  starting  values  in  a  program  that  varies  the  crystal-field 
parameters  to  obtain  a  best  least-squares  fit  to  the  experimental  data.  The 
resulting  crystal-field  parameters  are  used  in  conjunction  with  odd-fold 
lattice  sums  to  obtain  the  Judd-Ofelt  intensity  parameters,  the  level-to-level 
line  strengths  for  electric-  and  magnetic-dipole  transitions,  and  the  branch¬ 
ing  ratios  for  the  entire  rare-earth  series.  Comparisons  to  corresponding 
calculations  on  other  host  materials  are  made  throughout. 


*M.  J.  Weber,  Handbook  of  Lasers,  Chemical  Rubber  Company,  Cleveland.  Ohio 
(1971),  chapter  13,  p  371. 

2C.  A.  Morrison  and  R.  P.  Leavitt,  Spectroscopic  Properties  of  Triply 
Ionized  Lanthanides  in  Transparent  Host  Materials,  in  Handbook  on  the  Physics 
and  Chemistry  of  Rare  Earths ,  Vol  5,  K.  A.  Gschneidner  and  L.  Eyring,  Jr., 
eds.,  North  Holland,  New  York  (1982),  p  632. 

3A.  A.  Kaminskii ,  Kh.  S.  Bagdasarov,  G.  A.  Bogomolova,  M.  M.  Gritsenko,  A. 
M.  Kevorkov ,  and  S.  E.  Sarkisov ,  Luminescence  and  Stimulated  Emission  of  Nd^+ 
Ions  in  Gd ^Sc 2<3a ^0 ^  Crystals,  Phys.  Stat.  Sol.  (a)  34  (1976),  K109. 

4A.  A.  Kaminskii,  Laser  Crystals,  Springer-Verlag ,  New  York  (1981),  p  133. 
5D.  Pruss,  G.  Huber,  A.  Biemowski,  V.  V.  Laptev,  I.  A.  Shcherbakov,  and  Y. 
V.  Zharikov,  Efficient  Cr3  Sensitized  Nd3  :GdScGa — Garnet  Laser  at  1.06  ym, 
Appl .  Phys.  B28  (1982),  358.  See  also  E.  V.  Zharikov,  V.  V.  Laptev,  E.  I. 

Sidorova,  Yu.  P.  Timefeev,  and  I.  A.  Shcherbakov,  Absolute  Quantum  Efficiency 
of  the  Luminescence  of  Cr3  Ions  in  Gadolinium  Gallium  and  Gadolinium  Scandium 
Gallium  Garnet  Crystals,  Sov.  J.  Quant.  Electron.  12_  (1982),  1124. 

Dutoit ,  J.  C.  Walling,  and  D.  S.  Boudreaux,  Microwave  Phonon  Attenuation 
in  Rare-Earth  Garnets:  Ion-Phonon  Interactions,  Phys.  Rev.  Bll  ( 1975) ,  34 9. 

7S.  H.  Wemple  and  W.  J.  Tabor,  Refractive  Index  Behavior  of  Garnets,  J. 
Appl.  Phys.  44_  (1973),  1395. 

8C.  A.  Morrison  and  R.  P.  Leavitt,  Crystal  Field  Analysis  of  Triply  Ionized 
Rare  Earth  Ions  in  Lanthanum  Trifluoride,  J.  Chem.  Phys.  7_1_  (1979),  2366. 

9R.  P.  Leavitt  and  C.  A.  Morrison,  Crystal  Field  Analysis  of  Triply  Ionized 
Rare  Earth  Ions  in  Lanthanum  Trifluoride.  II.  Intensity  Calculations,  J. 
Chem.  Phys.  73_  (1980),  749. 
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2.  CRYSTAL  FIELD 


The  crystallographic  data10  on  GdjSc2Ga;30.j  2  are  given  in  table  1.  We  were 

unable  to  find  a  refined  x-ray  study  of  this  material  and  have  taken  the 

oxygen  positions  to  be  the  same  as  reported  for  Y^Ga^O^*  The  lattice  sums, 

Anm,  were  evaluated  using  these  x-ray  data.  Effective  charges,  Z^  (q^  =  eZi ) , 

in  units  of  the  magnitude  of  the  electron  charge,  were  taken  as  Zq^  =  ZSc  =  3, 

=  1,  and  ZQ  =  -1.5.  The  resulting  point -charge  contributions  to  the 

A  are  shown  in  column  2  of  table  2.  The  point-dipole11  and  self-induced1^ 
nm  , 

contributions  to  the  A  were  also  evaluated  and  are  shown  in  columns  3  and 

nm 

4.  In  the  evaluation  of  these  latter  contributions,  the  polarizability  of 
oxygen  was  taken  to  be  one -tenth  of  the  handbook  value  (based  on  previous 
experience).  The  point-charge  contributions  to  the  Anm  for  the  isostructural 
materials  Y^Al^O^13  and  Y^Ga^O^14  are  included  in  table  3  for  comparison. 


TABLE  1.  CRYSTALLOGRAPHIC  DATA  FOR  GdjScjGajO^ 
(CUBIC  SPACE  GROUP  Ia3d,  No.  230,  Z  =  8)a 


Ion 

Pos ition 

Symmetry 

X 

y 

z 

Gd 

24(c) 

D2 

0 

1/4 

1/8 

Sc 

16(a  ) 

C3i 

0 

0 

0 

Ga 

24(d) 

S4 

3/8 

0 

1/4 

0 

96(h) 

C1 

-0.0272 

0.0558 

0.1501 

aFractional  oxygen  positions  are  from  isostructural  material  Y^Ga^O^ 
(Morrison  and  Leavitt,  ref  2,  p  643).  Lattice  constant  is  a  =  12.5668 

(Brandle  and  Barns,  ref  10). 


10C.  D.  Brandle  and  R.  L.  Barns,  Crystal  Stoichiometry  and  Growth  of  Rare- 
Earth  Garnets  Containing  Scandium,  J.  Cryst.  Growth  20_  (1979),  1. 

1JC.  A.  Morrison,  Dipolar  Contributions  to  the  Crystal  Field  in  Ionic 
Solids,  Solid  State  Commun.  18_  (1976),  153, 

l2C.  A.  Morrison ,  G,  F.  de  S<3,  and  R.  P,  Leavitt,  Self -Induced  Multipole 
Contribution  to  the  Single-Electron  Crystal  Field,  J,  Chem.  Phys.  76_  (1982), 
3899. 

E.  Wortman,  C.  A.  Morrison,  and  N.  Karayianis,  Rare  Earth  Ion-Host 
Lattice  Interactions:  11.  Lanthanides  in  YjAlrO.j,  Harry  Diamond 

Laboratories ,  HDL-TR-1773  (August  1976). 

lkN.  Karayianis,  D.  E.  Wortman,  and  C.  A.  Morrison,  Rare  Earth  Ion-Host 
Lattice  Interactions :  12.  Lanthanides  in  YjGas Gi2'  Harry  Diamond 

Laboratories,  HDL-TR-1793  (July  1977). 
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where  the  pn(Nd)  are  radial  factors  given  in  an  earlier  work.15  These  Bnm 
were  then  used  as  initial  values  in  a  least-squares  fit  to  the  experimental 
data  on  NDiGdjScjGa-jC^  2«  The  best  fit  for  all  the  reported  data  was  extremely 
poor,  and  the  resultant  Bnm  were  far  from  those  obtained  for  Nd3  +  in  the 
isostructural  compounds  Y3A15°i2  antl  Y3Ga5°i2*  Tllese  results  are  shown  in 


TABLE  4.  BEST  FIT  OF  CALCULATED  ENERGY  LEVELS  TO  EXPERIMENTAL 
ENERGY  LEVELS  OF  Nd : Gd jSc 2Ga 30, 2  USING  ALL  EXPERIMENTAL  DATA 
OF  KAMINSKI  I  ET  AL  (ref  3,  4)  (E  =  experimental,  T  *  theoretical)3 


Term 

Prefix1* 

Energy 

(cm-1  ) 

9/2 

E 

- 

0 

107 

168 

263 

763 

T 

* 

-84 

136 

201 

304 

746 

*1/2 

E 

1 

978 

010* 

069 

109 

393 

431 

T 

2 

000 

023 

067 

M2 

339 

443 

^x 

^  1  3/2 

E 

3 

907 

920 

000* 

010 

025 

380 

421 

T 

3 

890 

900 

933 

013* 

192 

337 

406 

4  x 

*15/2 

E 

5 

111 

812 

914 

959 

494* 

510 

557 

T 

5 

772 

812 

915 

009* 

431 

474 

591 

E 

6 

647 

T 

6 

664 

4p 

3/2 

E 

T 

11 

1 1 

434 

4  34 

499 

494 

aRms  deviation 

is  45 

.  95  cm- 1 . 

Resulting  B  _  fern-1 

)  are  B20  = 

187 ,  &22  = 

-773,  B4(j  = 

-390, 

®42  = 

-475,  B4i} 

=*  -21 57, 

B60 

=  -2081 

'  b62  = 

-307 , 

B64  =  706, 

and  Bgg  -  -211. 

“Prefix  for  entries  to  right. 

Asterisk 

indicates  that,  for  that 

ent  ry  and 

succeeding 

entries  in 

that  row 

,  prefix 

is 

increased  by 

one. 

Thus  for 

experimental  data 

in  ‘*1 

11/2  we  ,uve 

I  975  010*  0 69  ...  =  1978  2010  20 69  ... 


15C.  A.  Morrison,  N.  Karayianis,  and  D.  E.  Wortman,  Rare  Earth  Ion-Host 
Lattice  Interactions:  4.  Predicting  Spectra  and  Intensities  of  Lanthanides  in 
Crystals,  Harry  Diamond  Laboratories,  HDL-TR-1816  (June  1977). 
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For  several  reasons,  we  then  chose  to  remove  the  level  at  4000  cm-*  in  an 

attempt  to  remedy  the  poor  fit.  First,  our  lattice-sum  calculation  predicts 

no  line  in  the  vicinity  of  the  4000-cm-*  line.  Second,  examination  of  the 

data  on  Nd3+  in  other  garnets2 3  reveals  at  most  four  lines  below  4200  cm-* 

within  the  4 1 13/2  manifold.  Third,  Kaminskii  et  al^  mention  extra  lines  in 

their  data.  Finally,  and  most  importantly,  rejection  of  the  4000-cm-*  line 

leads  to  a  vastly  improved  fit  of  the  data.  The  best  fit  to  the  experimental 

data  (minus  the  4000-cm-*  line)  was  obtained  again;  results  of  this  fit  are 

shown  in  table  5.  The  rms  deviation  between  the  calculated  and  experimental 

energy  levels  was  very  small  (~2  cm-* )  and  represents  one  of  the  best  fits 

that  we  have  ever  obtained  to  experimental  energy-level  data  in  crystals. 

Shown  in  table  6  are  the  best-fit  Bnm  and  the  phenomenological  Anm  (for  even 

n)  obtained  from  the  best  fit  B„  by 

nm  J 

Anm  =  Bnn/Pn<Nd>  '  <2> 

where  pn(Nd)  is  the  same  as  in  equation  (1). 


TABLE  5.  BEST  FIT  OF  CALCULATED  ENERGY  LEVELS  TO  EXPERIMENTAL 
ENERGY  LEVELS  OF  Nd:  GdjSCjGa-jO,  2  WITH  LEVEL  AT  4000  cm"1 
REMOVED  FROM  EXPERIMENTAL  DATA  GIVEN  IN  TABLE  4a 


Term 

Prefi 

xb 

Energy  (cm 

') 

4  j 

A9/2 

E 

- 

0 

107 

168 

263 

763 

T 

- 

-1 

105 

168 

267 

762 

4t 

A1  1  /2 

E 

1 

978 

010* 

069 

109 

39  3 

431 

T 

1 

979 

006* 

070 

107 

394 

4  35 

4. 

1 1  3/2 

E 

3 

907 

920 

010* 

025 

_ 

380 

421 

T 

3 

911 

923 

009* 

020 

367 

381 

419 

4t 

A  1  5/2 

E 

5 

777 

812 

914 

959 

494* 

510 

557 

T 

777 

809 

916 

962 

49S* 

S 1  2 

557 

E 

6 

647 

T 

6 

643 

4F 

r  3/2 

E 

1 1 

434 

499 

T 

1 1 

435 

498 

aRms  is 

2.37  cm~ 

l  . 

Resulting 

best -fit 

Brm  *re 

given 

in  tab! 

e  6 . 

For  explanation 

of 

notation , 

see  footnotes  of 

table 

4. 

2C.  A.  Morrison  and  R.  P.  Leavitt,  Spectroscopic  Properties  of  Triply 
Ionized  Lanthanides  in  Transparent  Host  Materials,  in  Handbook  on  the  Physics 
and  Chemistry  of  Rare  Earths,  Vol  5,  K.  A.  Gschneidne r  and  L.  Eyring,  Jr., 
eds.,  North  Holland,  New  York  (1982),  p  6 32. 

3A.  A.  Kaminskii,  Kh.  S.  Bagdasarov,  G.  A.  Bogomolova,  M.  M.  Gritsenko,  A. 
M.  Kevorkov,  and  S.  E.  Sarkisov,  Luminescence  and  Stimulated  Emission  of  JVd3+ 

Ions  in  Gd3Sc2Ga30j2  Crystals,  Phys.  Stat.  Sol.  (a )  34  (1976),  K109. 


TABLE  6.  BEST  FIT  B  FOR  Nd  IN  Gd ,Sc,Ga ->0. ,  AND 

flirt  -3  4  i  '  4 

PHENOMENOLOGICAL  Anm  (Anm  =  Bnn/Pn<Nd)) 


nm 

B  (cm  1  ) 
nm 

A  (cm  V) 
nm 

20 

434 

2544 

22 

90 

529 

40 

-67 

-115 

42 

-1818 

-3148 

44 

-921 

-1595 

60 

-1485 

-934 

62 

-679 

-427 

64 

771 

485 

66 

-661 

-416 

The  Anm  obtained  by  using  equation  (2),  as  given  in  table  6,  were  then 

compared  to  the  various  contributions  to  the  Anm  given  in  table  2.  It  is 

apparent  that  if  different  polarizabilities  and  effective  charges  were  chosen, 

the  values  of  Anm  of  table  2  could  be  brought  into  better  agreement  with  those 

of  table  6.  However,  because  of  the  lack  of  refined  x-ray  data  for 

GdjSc 2Ga3012  (that  is,  the  oxygen  positions  are  unknown),  no  attempt  was  made 

to  investigate  such  changes.  The  phenomenological  Anm  of  table  6  for  even  n 

were  used  in  equation  (1)  to  obtain  Bnm  for  the  entire  rare-earth  series. 

These  B  are  given  in  table  7. 
nm  ^ 


TABLE  7.  SMOOTHED  CRYSTAL-Fl  ELD  PARAMETERS  B  (cm'1)  FOR  ALL 

-5 .  nm 

RARE-EARTH  IONS  R FOR  R : Gd jSc 2Ga 30 1 2 


R 

f,2i1 

B22 

B40 

B42 

«44 

0 

v£ 

E 

*62 

B64 

B66 

Ce 

46ft 

97.  S 

-86.  0 

-2372 

-1202 

-2187 

-1000 

1  1  36 

-974 

Pr 

447 

9.3.0 

-74.  4 

-2035 

-1031 

-1  752 

-801 

910 

-780 

N  1 

414 

90.  3 

-66.5 

-1818 

-921 

-1485 

-679 

771 

-661 

Phi 

42  7 

ftft .  9 

-61.5 

-1680 

-851 

-1  328 

-607 

690 

-591 

■4n 

4  24 

ftft.  3 

-SR.  1 

-1589 

-805 

-1234 

-564 

641 

-549 

Eu 

424 

ftfl .  2 

-55.  7 

-1522 

-771 

-1  168 

-534 

606 

-520 

424 

Rft.  3 

-53.6 

-1465 

-742 

-1109 

-507 

576 

-494 

Th 

426 

ftft.  6 

-51 . 7 

-141  3 

-716 

-1049 

-480 

545 

-467 

riy 

42R 

89.  0 

-50.0 

-1366 

-692 

-991 

-453 

515 

-441 

Ho 

410 

89.  6 

-48.  6 

-1  327 

-672 

-945 

-432 

491 

-421 

Er 

4  34 

90.  3 

-47.  5 

-1299 

-658 

-918 

-420 

477 

-409 

Tm 

43ft 

91 . 2 

-46.  7 

-1276 

-646 

-901 

-41  2 

468 

-401 

Yh 

442 

91 . 9 

-45.  3 

-1239 

-628 

-852 

-390 

442 

-3  79 

The  odd- n  A„  given  in  column  2  of  table  2  were  used  to  calculate  the 
nm 

Judd-Ofelt  intensity  parameters,  for  all  the  rare-earth  ions  in 

Gd-jSc2Ga301 2>  The  are  given  in  table  8.  In  table  9  the  corresponding  (2^ 
are  calculated  for  R:Y3A15012  and  in  table  10  the  are  given  for  RtLiYF^. 
The  results  in  tables  9  and  10  are  included  here  for  comparison. 
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TABLE  8.  CALCULATED  JUDD-OFELT  INTENSITY  PARAMETERS 
FOR  R:Gd3Sc2Ga3012d 


fi2  (10 


Q  (10'20cm2) 

4 


-20  2 
ii6(10  cm  ) 


Ce 

12.75 

2.81  7 

26.  21 

Pr 

7.136 

1.423 

11.19 

Nd 

6.880 

0.8867 

6.  1  39 

Pm 

3.982 

0.6923 

4.  717 

Sm 

3.472 

0.5808 

3.903 

Eu 

2.  747 

0.4554 

2.  764 

Gd 

2.150 

0. 3540 

1.909 

Tb 

3.  754 

0.5356 

4. 615 

Dy 

2.761 

0. 3963 

2.866 

Ho 

2.  255 

0. 3210 

2.093 

Er 

2.  162 

0. 2998 

1.974 

Tm 

2.103 

0. 2846 

1. 906 

Yb 

1.774 

0.  2360 

1.452 

a0dd~n 

A_,_  used  in  calculation  of  are  from  column  2  of 

nm  k 

table  2. 

TABLE  9.  CALCULATED  JUUD-OFEl.f  INTENSITY 

PARAMETERS 

FOR  R  :  Y-jAl^Oj  .,■* 

R 

-20  2 
*2(10  cm  ] 

i  si  .(  10"2ncmJ 

4 

> 

-20  2 
.  :fe(  10  cm  ) 

Ce 

1.933 

8.  751 

39.  78 

Pr 

1.082 

4.413 

16.96 

Nd 

1.043 

2.855 

9.  293 

Pm 

0.60  36 

2.  197 

7.  147 

Sm 

0.5264 

1 . 871 

5.864 

Eu 

0.4165 

1 . 422 

4.  186 

Gd 

0. 3259 

1.065 

2.  887 

Tb 

0.5692 

2.066 

7.029 

Dy 

0.4186 

1 . 422 

4.  360 

Ho 

0. 3418 

1.112 

3.  181 

Er 

0.  3278 

1.057 

3.001 

Tm 

0.  3188 

1.023 

2.  899 

Yb 

0. 2690 

0. 8288 

2.  208 

,3*'>dd-foZd  A  used  in  calcu 

nn 

4  df  ion  of  ^  are  from  c- 

C  u. fi;i  ?.  of 

tiblo  ?. 

odd  A 

nm 

of  column  2  of 

table  2  were  also 

used  to 

calculate 

to-line  transition  probabilities  for  electric-dipole  and  magnetic-dipole 
transitions  for  all  the  rare-earth  ions,  R3+,  in  R: Gd3Sc2Ga301 2.  Also,  the 
branching  ratios  and  g-factors  were  calculated.  These  results  are  included  in 
the  form  of  microfiche  copies  of  computer  outputs,  at  the  end  of  the  report; 
these  computer  outputs  are  explained  in  appendix  A.  The  branching  ratios 
for  *4f’3/2  to  ^Ij  for  J  =  15/2,  13/2,  11/2,  and  9/2  are  shown  in  figure  1 
for  Nd3+  in  Y3A15°i2  <YAG)«  LiYF^{  YLF ) ,  and  GdjSCjGa^  2.  Results  for  Nd3+  in 
YAG  and  YLF  are  included  because  experimental  data  are  available  for  these 
host  materials.  It  is  apparent  that  the  calculated  results  agree  well  with 
the  experimental  results;  this  gives  further  confidence  in  the  accuracy  of  the 
calculations  on  Gd3Sc2Ga303 2<  Thus,  if  the  x-ray  positions  assumed  for 
Gd3Sc2Ga3012  are  a  good  approximation,  we  should  expect  the  results  shown  in 
figure  1  for  Nd3+  in  this  host  to  be  quite  accurate. 


TABLE  10.  CALCULATED  JUDD-OFELT  INTENSITY  PARAMETERS 
FOR  R:LiYF4a 


R 

- TTfi - 5 - 

ii2(10  cn»  ) 

- — 5 - 

1K<10  *  cm  ) 

4 

-20  2 
tf6<10  cm  ) 

Ce 

22.28 

6. 161 

52.  72 

Pr 

12.47 

3.111 

22.48 

Nd 

12.02 

1.949 

12.  33 

Pm 

6.956 

1.519 

9.476 

Sm 

6.066 

1.277 

7.845 

Eu 

4.  799 

0.  9967 

5.551 

Gd 

3.755 

0.7706 

3.830 

Tb 

6.559 

1.213 

9.  304 

Dy 

4.824 

0.6864 

5.773 

Ho 

3.939 

0.7139 

4.213 

Er 

3.777 

0.6686 

3.975 

Tm 

3.673 

0.6367 

3.839 

Yb 

3.  100 

0.5260 

2.924 

aEven-n  B  are  from  D.  E.  Uortman,  N.  Karay ianis ,  and  C.  A.  Morrison, 
Rare  Earth  Ion-Host  Lattice  Interactons:  6.  Lanthanides  in  LiKF 4,  Harry 
Diamond  Laboratories,  HDL-TR-1770  (August  1976),  Odd-n  Anff(  are  from  L. 
Esterowitz,  F.  J,  Bartoli,  R,  E.  Allen,  D.  E.  Uortman,  C.  A.  Morrison, 
and  R.  P.  Leavitt,  Energy  Levels  and  Line  Intensities  of  Pr3  in  LiYF 4, 
Phys.  Rev.  B19  (1979),  6442. 


W„  uyf4 


Figure  1 .  Multiplet-to-multiplet  fluorescence  branching  ratios  (in  percent) 
for  transitions  originating  at  the  ^2/2  ^eve^  f°r  in  various 

hosts: 

a,  calculated,  K.  K.  Deb,  R.  G.  Buser,  C.  A.  Morrison,  and  R.  P.  Leavitt,  Crystal  Fields  and 
Intensities  of  Triply  Ionized  Rare-Earth  Ions  in  Cubic  Lanthanum  Oxyfluoride:  An  Efficient  4Fj.- 
♦  4I9/2  LaOF : Nd  Laser,  J.  Opt.  Soc.  Am.  71_  (1981),  1463; 

b,  experimental,  R.  K.  Matts,  Jr.,  Branching  Ratios  for  YAlG:Nd3*,  J.  Opt  Soc.  Am.  ^  (1971),  123; 

c,  experimental,  S.  Singh,  R.  G.  Smith,  and  L.  G.  Van  Uitert,  Stimulated  Emission  Cross  Section  and 
Fluorescent  $iantum  Efficiency  of  Nd3*  in  Yttrium  Garnet  at  Room  Temperature,  Phys.  Rev.  BIO 
(1974),  2566; 

d,  experimental,  H.  P.  Jenssen,  Center  for  Materials  Science  and  Engineering  and  Department  of 
Electrical  Engineering,  Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139  (personal 
communication ) ; 

e,  experimental,  T.  L.  Lomheim  and  L.  G.  DeShazer,  New  Procedure  of  Determining  Neodymium 
Fluorescence  Branching  Ratios  as  Applied  to  25  Crystals  and  Glass  Hosts,  Opt.  Commun.  _24  (1978), 
89; 

f,  calculated,  this  work. 


3.  SCANDIUM  AND  GALLIUM  SITES 


To  investigate  the  possibility  of  other  site  occupancy  by  Nd3  +  or  other 
ions  in  Gd3Sc2Ga30,j the  distances  between  ions  and  the  number  of  ion  pairs 
at  a  given  distance  were  calculated  for  all  the  sites  in  this  compound.  The 
results  are  given  in  table  11.  The  ionic  radii  of  several  ions  in  different 
coordination  are  given  in  table  12.  The  coordination  number  (CN),  the  number 
of  nearest  ions,  for  the  Gd  site  is  assumed  to  be  8,  since  the  two  groups  of 
oxygen  are  nearly  the  same  (see  table  11)  and  the  next  set  of  ions  is  much 
further  away  (~3.14  A).  The  ionic  radii  of  Nd3+  and  Gd3+  for  CN  =  8  are 
nearly  the  same,  so  we  should  expect  the  Nd3+  to  predominantly  enter  the  Gd3+ 
site  in  Gd^ScjGajO.!  However,  it  is  possible  that  some  of  the  Nd3+  may  enter 
the  Sc3  +  site  in  this  compound.  For  this  reason,  the  point-charge  lattice 
sums  Anm  were  calculated  for  the  Sc3+  site  in  Gd3Sc2Ga301 2;  the  results  are 
given  in  table  13.  Also  included  in  table  13  are  the  Bnm  for  Nd3+  calculated 
using  equation  (1).  In  these  calculations,  the  effective  charges  on  the  ions 
are  the  same  as  used  in  the  calculation  of  table  2.  The  crystal-field 
parameters,  Bnm>  for  Nd3+  in  the  scandium  site  as  given  in  table  13  are  un¬ 
realistically  large.  This  possibly  is  due  to  the  inaccuracies  in  the  x-ray 
data.  No  further  calculation  was  done  using  the  Bnm  for  the  scandium  site. 


TABLE  11.  DISTANCE  TO  NEAREST  OXYGEN  IONS, 

R  (A),  AND  NUMBER  OF  OXYGEN  IONS  AT  THAT  DISTANCE  FOR  METAL 
SITES  IN  Gri  3Sc2Ga  jO,  2a 


Site 

R 

Number 

R 

Number 

Gd 

2.  30 

4 

2.48 

4 

Sc 

2.04 

f> 

- 

- 

Ga 

1.89 

4 

- 

- 

aX-ra y  data  used  are  given  in  table  l. 


TABLE  12.  IONIC  RADII  (A)  AND  COORDINATION  NUMBERS  (CN )  OF 
SEVERAL  IONS'1 


Ion 

Ionic  ra»1ius 

CN 

Gd3+ 

vC 

’ 

8 

Nd  3+ 

1.12 

8 

He3* 

.  87 

6 

0a3+ 

.47 

4 

Cr  3+ 

.615 

3R.  D.  Shannon  and  C.  .7.  Prewitt,  Effective  Ionic  Radii  in  Oxides  and 
Fluorides,  Acta .  Crys  t.  R2S  (1060),  121.  Ionic  radius  is  based  on 

oxygen  ion  in  6  coordination  having  value  1.40  A. 

°Ionic  radius  of  Cr  *  in  4  coordination  is  not  given  in  a  hive 
reference. 
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TABLE  13.  POZ  NT -CHARGE  LATTICE  SUMS  An>  (ciT’/A  >  FOR 
Sc3*  SITE  IN  Gd 3Sc  jGa 30 ,  2 


B  (Nd)b 

nm 


R 

63 

2099 

33 

i 

63 

52 

i 

R 

66 

687 

10' 

Z 

66 

25 

’’x-ray  data  used  are  given  in  table  1.  Effective  charges  on  ions  are 
as  in  table  2.  Lattice  sum  was  performed  for  Sc3*  site  at  (0,0,0),  and 
a  rotation  was  performed  so  that  Cjj  axis  is  along  a  body  diagonal, 
Siller  angles  used  are  a  ”  -45*,  0  *  -54.736“,  and  y  •  28.319°. 


bB^lNd)  are  calculated  using  agnation  ( 1). 


The  point-charge  lattice  sums  were  calculated  for  the  Ga  site  (S^ 
symmetry)  in  Gd3Sc2Ga301 2»  the  results  are  given  in  table  14.  It  is 
highly  unlikely  that  the  Nd3  +  enters  this  site,  because  of  the  large  mismatch 
in  ionic  radii  as  given  in  table  12.  However  it  is  quite  possible  that  Cr3  + 
doped  in  this  crystal,  as  reported  elsewhere,5  may  enter  this  site  or  perhaps 
the  Sc3+  site  (see  table  12).  If  the  Cr3  +  were  to  enter  the  Ga  site  substitu- 
tionally,  the  absorption  should  be  quite  intense  because  of  the  odd-n 
present  in  S4  symmetry.  Thus  it  would  be  possible  to  easily  pump  the  Cr3  ion 
and  observe  the  transfer  from  Cr3+  to  a  substitutional  rare-earth  impurity  in 
Gd3Sc2Ga301 2.  In  future  work,  we  expect  to  use  the  results  in  table  14  to 
estimate  the  splittings  of  levels  of  Cr3+  in  this  host. 


TABLE  14.  POINT -CHARGE  LATTICE  SUMS 

A_  (cm_,/An>  FOR  Ga3*  SITE  IN 
n  in 

Gd3Sc2Ga3°12 


aX-ray  data  used  are  given  in  table 
1.  Effective  charges  on  ions  are  as  in 
table  2.  Lattice  sum  was  performed  for 
Ga 3*  site  at  (0,0,0),  and  a  rotation  was 
performed  so  that  S ^  axis  is  along  z 
axis.  Euler  angles  used  are  a  •  0° ,  0  * 
90°,  and  y  -  59.5626“. 


5D.  Pruss,  G.  Huber,  A.  Biemowski ,  V.  V.  Laptev,  I.  A.  Shcherbakov,  and  Y . 
V.  Zharikov,  Efficient  Cr3+  Sensitized  Nd*+  :GdScGa — Garnet  Laser  at  1.06  uni/ 
Appl.  Phys.  B28  (1982),  358.  See  also  E.  V.  Zharikov,  V.  V.  Laptev,  E.  I. 
Sidorova,  Yu.  P.  Timefeev,  and  I.  A.  Shcherbakov,  Absolute  Quantum  Efficiency 
of  the  Luminescence  of  Cr3+  Ions  in  Gadolinium  Gallium  and  Gadolinium  Scandium 
Gallium  Garnet  Crystals,  Sov.  J.  Quant.  Electron.  12  (1982),  1124. 


4.  CONCLUSION 


We  have  analyzed  the  experimental  spectra  of  Nd3  +  in  the  host 
GdjSCjGa^Oj  2*  By  omitting  one  of  the  reported  energy  levels  we  were  able  to 
obtain  an  excellent  fit  (rms  ~  2  cm-M  to  the  experimental  data.  The  re¬ 
sulting  phenomenological  crystal-field  parameters,  Bnm,  were  used  to  obtain 
Bnm  for  the  entire  rare-earth  series.8  Point-charge  lattice  sums  Anm  were 
calculated  for  effective  charges  of  ZQ  =  -1.5,  ZGa  =  +1,  Zgc  =  +3,  and  ZM  = 
+3.  The  even-n  Bnm  calculated  as  in  earlier  work2  agreed  quite  reasonably 
with  the  phenomenological  Bnm  obtained  by  fitting  the  experimental  data.  The 
odd-n  Anm  from  the  point-charge  lattice  sums  were  used  along  with  the  phenome¬ 
nological  Bnn)  to  calculate  the  crystal -field-split  energy  levels,  transition 
intensities,  branching  ratios,  g  values,  and  lifetimes  for  the  low-lying 
energy  levels  of  all  the  rare-earth  ions  in  the  Gd3  +  sites  in  Gd3Sc2Ga-j01 2- 
The  resulting  branching  ratios  for  Nd3+  and  several  other  computed  quantities 
were  compared  with  standard  laser-host  materials  (YAG,  YLF)  with  few  striking 
differences.  The  laser  advantages  of  the  garnet  Gd3Sc2Ga3012  appear  to  lie  in 
the  fact  that  the  material  can  be  easily  doped  with  transition  elements  such 
as  Cr3+  and  still  produce  high-optical-quality  crystals. 
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APPENDIX  A 


APPENDIX  A. — COMPUTER  OUTPUT 

Detailed  intensity  ( line-to-line )  and  g-f actor  calculations  have  been 
performed,  and  the  results  are  attached  in  the  form  of  microfiche  copies  of 
computer  printouts  for  each  ion  in  the  lanthanide  series,  from  Ce3+  to  Yb3+. 
In  this  appendix,  we  describe  the  contents  of  the  computer  printouts  and 
present  excerpts  from  them. 

The  computer  printout  for  each  ion  presents  the  following  information,  in 
sequence,  as  follows: 

(1)  even-k  crystal-field  parameters,  centroid  parameters,  and  crystal- 
field-split  energy  levels; 

(2)  approximate  Judd-Ofelt  intensity  calculations  (no  J  mixing)  including 
the  Judd-Ofelt  parameters,  multiplet-to-multiplet  line  strengths,  lifetimes, 
and  branching  ratios ; 

(3)  full  Judd-Ofelt  squared  matrix  elements  for  u  (E  parallel  to  D2  axis) 
and  a  (E  perpendicular  to  D2  axis)  between  all  crystal-field-split  levels; 

(4)  a  summary  of  the  full  Judd-Ofelt  calculations  (J  mixing  included), 

including  multiplet-to-multiplet  line  strengths,  lifetimes,  and  branching 

ratios ; 

(5)  approximate  magnetic -dipole  intensity  calculations  (no  J  mixing), 

including  multiplet-to-multiplet  line  strengths,  lifetimes,  and  branching 

ratios ; 

(6)  full  magnetic -dipole  squared  matrix  elements  for  a  (B  parallel  to  D2 
axis)  and  n  (B  perpendicular  to  D2  axis)  between  all  crystal-field-split 
levels; 

(7)  a  summary  of  the  full  magnetic-dipole  calculations  (J  mixing 
included),  including  multiplet-to-multiplet  line  strengths,  lifetimes,  and 
branching  ratios ;  and 

(8)  g  factors  (for  an  odd  number  of  f  electrons). 

An  excerpt  of  the  energy  level  table  for  Pr3  +  is  given  in  table  A-1 .  The 
crystal-field-split  levels  are  numbered  in  sequence,  beginning  with  the  ground 
state  (these  sequence  numbers  are  used  later  as  identifiers  in  the  line-to- 
line  intensity  calculations).  The  "PCT.  PURE”  column  gives  an  indication  of 
the  amount  of  J  mixing  by  the  crystal  field.  The  crystal  quantum  number 
column  ("2MU")  identifies,  for  an  even  number  of  f  electrons,  the  irreducible 
representation  of  C2  to  which  a  level  belongs;  levels  with  2u  =  0  are  r1  and 
levels  with  2jj  =  2  are  r 2*  Eor  odd  number  of  electrons,  all  levels  are 
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Kramers  doublets  with  2 vi  =  1  (r3  4).  The  theoretical  energy  levels  ("THEOR. 
ENERGY")  correspond  to  those  obtained  using  aqueous  free-ion  parameters.1 
(The  "EXP.  ENERGY"  column  is  not  relevant  to  this  discussion.  ) 


TABLE  A-1  .  SAMPLE  OF  ENERGY  LEVELS  TABLE  FROM 
COMPUTER  OUTPUT  FOR  Pr 3+  IN  Gd  jSCjGa  jO,  2 


FREE  ION 

PCT 

PURE 

2MU 

THEOR. 

ENERGY 

EXP. 

ENERGY 

1  3H  4 

99.9 

0 

0.0 

0.0 

2  3H  4 

99.8 

0 

38.4 

0.0 

3  3H  4 

99.6 

2 

44.8 

0.0 

4  3H  4 

98.5 

0 

452.8 

0.0 

5  3H  4 

96.1 

0 

474.6 

0.0 

6  3H  4 

95.8 

2 

507.4 

0.0 

7  3H  4 

95.7 

2 

516.4 

0.0 

8  3H  4 

94.4 

2 

717.2 

0.0 

9  3H  4 

94.6 

0 

731 .2 

0.0 

64  ID  2 

99.9 

2 

16733.5 

0.0 

65  ID  2 

99.9 

0 

16746. 3 

0.0 

66  ID  2 

99.9 

0 

17273.2 

0.0 

67  ID  2 

99.9 

0 

17438.5 

0.0 

68  ID  2 

99.8 

2 

17458.4 

0.0 

A  sample  of  the  computer  printout  for  the  approximate  Judd-Ofelt 
calculations  for  Pr3  +  is  given  in  table  A-2.  In  the  table,  the  upper  right 
elements  give  the  multiplet-to-multiplet  line  strengths  in  units  of  square 
angstroms.  Thus,  for  instance,  the  total  line  strength  for  the  3H^  •»  3F-j 
transition  is  5.0  x  10'5  A2.  The  lower  left  elements  in  the  table  give  the 
branching  ratios  for  spontaneous  emission,  assuming  that  all  levels  within  a 
multiplet  are  equally  populated.  The  branching  ratios,  for  example,  for 
emission  from  the  ^2  level  are  as  follows: 


0.06 

for 

lD2 

*► 

3fV 

0. 0069 

for 

lD2 

-► 

3fv 

0.049 

for 

lD2 

3p 

2' 

0.081 

for 

lD2 

-► 

3V 

0.0063 

for 

1 D2 

-► 

3H5' 

0.72 

for 

lD0 

* 

3ha- 

For  a  given  initial  state,  all  the  branching  ratios  add  up  to  unity.  The 
calculation  of  the  branching  ratios  neglects  the  dependence  of  the  index  of 


W.  T.  Carnall,  P.  R .  Fields,  and  K.  Rajnak,  El  'ctronic  Levels  in  the 
Trivalent  Lanthanide  Aquo  Ions.  Pr^+,  Nd^+ ,  Sm3+,  Dy3*,  f/o3+,  Er3+,  and  Tm3+, 
J.  Chem.  Phys.  49^  (1968),  4412. 
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refraction  on  frequency.  The  diagonal  elements  in  the  table  give  the 
lifetimes  of  the  excited  multiplets,  again  assuming  equal  populations  of  the 
initial  sublevels.  The  quoted  lifetimes  are  in  seconds  and  should  be  divided 
by  n(n2  +  2)2/9,  where  n  is  the  index  of  refraction. 


TABLE  A-2.  SAMPLE  OP  COMPUTER  PRINTOUT  SUMMARY  OF  ELECTRIC- 
DIPOLE  TRANSITION  PROBABILITY  CALCULATIONS  FOR  Pr3+  IN 
GdjSc2Ga3012  USING  APPROXIMATE  JUDD-OFELT  THEORY 


3H  4 

3H  5 

3H  6 

IF  2 

3F  3 

3F  4 

1 G  4 

ID  2 

3H 

4 

0.0 

7.  IE-04 

1.6E-04 

1 . 9E-04 

8. 3E-04 

5.5E-04 

3.  IE-05 

6. IE-05 

3H 

5 

\ . oe+oo 

2. IE-01 

7. 5E-04 

7.8E-04 

5.0E-05 

5. 3E-04 

4.9E-04 

8.  IE-07 

3H 

6 

6. 3E-01 

3. 7E-01 

8. 9E-02 

3.  4E-04 

9. 9E-04 

6. 0E-04 

3.  2E-04 

1 . 7E-05 

3F 

2 

5.7E-01 

4. 3E-01 

2. 8E-03 

1 . 7E-02 

7.  3E-06 

1 . 0E-04 

9.  9E-06 

1 . 2E-05 

3F 

3 

9. 5E-01 

1.6E-02 

3.6E-02 

7. 3E-05 

4.4E-03 

1 .6E-05 

5.  9E-05 

2. 4E-06 

3F 

4 

7. 4E-01 

2. 2E-01 

3. 8E-02 

2. 3E-03 

6.4E-06 

5. 4E-03 

4.2E-04 

2. 4E-05 

1 G 

4 

9. 3E-02 

7. 0E-01 

1.6E-01 

3.4E-03 

7. 7E-03 

3. 6E-02 

4. IE-03 

9. 2E-05 

ID 

2 

7. 2E-01 

6. 3E-03 

8. IE-02 

4. 9E-02 

6.9E-03 

6.0E-02 

7.9E-02 

1 . 8E-03 

An  excerpt  from  the  squared  matrix-element  table  is  given  in  table  A-3  for 
the  x-spectrum,  •>  ^2  line-to-line  electric-dipole  transitions.  It  should 

be  kept  in  mind  that  the  polarizations  referred  to  in  this  and  the  following 
table  are  with  respect  to  the  individual  C2  axes  of  the  sites,  and  not  the  c- 
axis  of  the  crystal.  In  table  A-3,  one  number  is  presented  for  each 
transition,  the  squared  matrix  element  p2  of  the  z  component  of  the  electric 
dipole  operator,  in  units  of  square  angstroms. 


TABLE  A-3.  EXCERPTS  FROM  COMPUTER  OUTPUT  FOR  Pr 3+ : Gd jSc 2Ga jO, 2 
LINE-TO-LINE  71  ELECTRIC-DIPOLE  TRANSITION  PROBABILITY  TABLE,  FOR 
3H4  ►  ’d2  TRANSITIONS  WITH  2y  =  0  .  2y  =  0  AND  2y  »  2  ♦  2y  -  2 


2y  = 

o 

II 

£ 

♦ 

o 

65 

6667 

ID  2 

ID  21D  2 

1  3H 

4 

3. 568E-07 

5. 986E-06 

8.  124E-19 

2  3H 

4 

1 . 921E-18 

6. 741E-1 8 

6. 936E-08 

4  3H 

4 

8.  742E-19 

3.210E-20 

6. 083E-07 

5  3H 

4 

3.  344E-08 

5. 293E-07 

6.  175E-19 

9  3H 

4 

3. 1 36E-18 

1 . 374E-18 

4. 265E-06 

2y  *  2  ♦  2y  *  2 

64 

68 

ID  2 

in  2 

3  3H  4 

6.741E-06 

1  •  608:'-1  9 

6  3H  4 

4. 875E-07 

3.644E-18 

7  3H  4 

1.652E-17 

8.  572E-08 

8  3H  4 

2. 502E-21 

2.  1 34E-06 
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Table  A-4  shows  excerpts  from  the  squared  matrix  elements  table  for  the  o- 
spectrum,  ♦  1d2  electric-dipole  transitions.  Again,  the  polarization  (o) 

is  with  respect  to  the  site  D2  axis.  For  each  transition,  three  numbers  are 
presented  in  the  table:  p2y(max  >'  Pxy  (min),  and  0m. 

TABLE  A-4.  EXCERPTS  PROM  COMPUTER  OUTPUT  FOR  Pr  3+ :  Gd  jSc  jGa  }0,  2 
LINE-TO-LINE  0  ELECTRIC-DIPOLE  TRANSITION  PROBABILITY  TABLE,  FOR 
\  .  ’l)2  TRANSITIONS  WITH  2U  -  0  .  2y  =  2  AND  2y  =  2  ♦  2p  =  0 


2u  =  0  .  2g  =  2 


i 

2 

4 

5 

9 

3H  4 

3H  4 

3H  4 

3H  4 

3H  4 

f>4 

ID 

2 

4.686E-06 

2. 21 1 E-07 

3.  440E-09 

1. 146E-06 

6. 598E-08 

0.0 

0.0 

o 

o 

0.0 

o 

o 

i  68 

0.0 

0.0 

0.0 

0.  0 

o 

o 

ID 

2 

3.062E-06 

1 . 1 40E-07 

5. 466E-07 

2.  875E-07 

3.  126E-06 

4.  547E-I 3 

0.0 

5. 684E- 1 4 

2.842E-14 

o 

o 

90.0 

90.0 

o 

d 

90.0 

0.0 

>  2“  ■  2 

♦  2u  »  0 

65 

66 

67 

ID  2 

ID  2 

ID  2 

i 

1  3  3H  4 

2. 772E-07 

5. 163E-06 

3. 299E-06 

i 

2. 842E- 1 4 

4. 547E-1 3 

4.547E-13 

1 

1 

0.0 

0.0 

90.00 

6  3H  4 

t • 993E-09 

9.950E-08 

1 . 874E-07 

i 

1.11 0E-1 6 

0.0 

-2.842E-14 

0.0 

0.0 

90.00 

7  3H  4 

4.879E-07 

3.618E-06 

1 . 553E-06 

8.527E-14 

0.0 

0.0 

90.00 

90.00 

0.0 

8  3H  4 

2.  593E-06 

5.663E-07 

2. 482E-07 

0.0 

2.842E-14 

2.842E-14 

1 

90.00 

90.00 

0.  0 

The  remaining  computer  printout  for  intensity  calculations  consists  of 
multiplet-to-multiplet  tables  and  line-to-line  tables,  as  described  above. 
These  tables  have  the  same  format  as  those  already  described.  It  should  be 
remembered  that,  for  the  magnetic-dipole  tables,  it  and  o  are  reversed;  that 
is,  a  corresponds  to  the  B  vector  along  z,  and  it  corresponds  to  the  B  vector 
perpendicular  to  z.  Also,  for  the  magnetic-dipole  transitions,  the  radiative 
lifetime  is  to  be  divided  by  n^,  where  n  is  the  index  of  refraction. 

For  an  odd  number  of  electrons,  all  the  levels  are  Kramers  doublets,  and  g 
factors  have  been  computed.  The  g  tensor  can  be  characterized  by  four 
numbers:  g|  (B  parallel  to  z),  g^1  and  g^  (maximum  and  minimum  values 
of  g^  for  B  in  the  x-y  plane,  respectively),  and  9m  (the  angle  in  the  x-y 
plane  for  which  g^  is  a  maximum).  These  are  presented  on  the  last  page  of  the 
computer  output. 
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